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Abstract 

We consider lattice deformations (both continuous and topolog- 
ical) in the hexagonal lattice Hubbard model in the tight binding 
approximation to graphene, involving operators with the range up to 
next-to-neighbor. In the low energy limit, we find that these deforma- 
tions give rise to couplings of the electronic Dirac field to an external 
scalar (Yukawa) and gauge fields. The fields are expressed in terms of 
original defects. As a by-product we establish that the next-to-nearest 
order is the minimal range of deformations which produces the com- 
plete gauge and scalar fields. We consider an example of Stone-Wales 
defect, and find the associated gauge field. 



1 Introduction 

Graphene attracted a lot of interest in recent years from both theoretical and 
experimental communities due to its remarkable properties pQ. Especially 
intriguing feature of graphene's electronic structure is that in the low energy 
limit it is described by a (pseudo)relativistic massless Dirac fermion model [2]. 
This property provides a link between the condensed matter theory and high 
energy physics. Another interesting feature of graphene is its chirality. The 
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microscopic theory is non-relativistic and chiral non- invariant, but due to 
a phenomenon similar to fermion doubhng in lattice field theory, the chiral 
and Lorentz symmetries, as well as an internal symmetry emerge in the low 
energy (continuum) limit [3]. 

The nonabelian character of the internal symmetry of the system is in- 
spiring for a further quest into other links which graphene and high energy 
physics may have. In particular, an attractive idea is to discover the non- 
abelian gauge fields in graphene. A natural candidate for this role is provided 
by the phonon field or more generally by the defects of the hexagon lat- 
tice. Thus, experimentally it was discovered that the phonon couples to the 
electron wave function in an anomalous way. This lead to the suggestion, 
that this coupling should be realized, in fact, as a Dirac fermion/gauge or 
gravity coupling [SHH]. The phonon field can be regarded as a particular 
case of smooth defect. In more general case including topological defects, 
it was found that various types of lattice defects can be described in the 
low energy limit by coupling the fermionic field to various non-abelian field 
backgrounds P4TT] (see [12] for a review). 

The objective of this work is to consider most arbitrary deformation of the 
hexagonal Hubbard model, corresponding to the inclusion of various local or 
short range defects. The general problem related, in particular, to topological 
defects is that in many cases they do not admit, even in the low energy limit, 
a description in terms of continuous fields. In contrast, they may localize 
around several points of the momentum space. Our strategy in this work is 
to consider a generic setup and extract the modes coupled to the low energy 
electronic field. It appears, that it is the modes localized around the origin as 
well as around the Dirac points which are 'seen' by the low energy fermionic 
field. 

The plan of the remainder of the paper is as follows. In the next section 
we give a brief introduction to the low energy theory of the free electronic 
wave function and introduce the notations. In the following section we con- 
sider the deformation of the hexagonal lattice Hubbard model by inserting 
the terms into the Hamiltonian corresponding to local, nearest neighbor and 
next-to-nearest neighbor deformations of the model. We assume, that the 
deformation has an a priori arbitrary spectral composition. Then we take 
the low energy limit and select the modes which couple to the low energy 
fermionic wave function. The main result is summarized by eq. fl3.19p : the 
defects emerge as a U(2) scalar (Yukawa) and U(2) gauge couping. The re- 
sult also shows that the next-to-neighbor term is important to obtain the 
complete gauge/scalar coupling. At the same time this type of deformations 
is also saturating: adding more distant interaction terms to the Hamiltonian 
will not change the character of interactions in the low energy theory, but 
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just will renormalize it. We find these properties of the model remarkable. 
We also rule out the assumption existent in the literature that the coupling 
can be fermionic gravity type coupling. In the Section H] we consider an ap- 
plication of the general theory to the description of the Stone- Wales defect 
in terms of gauge fields. In addition we include two sections in the Appendix: 
One providing the complete information on the spinor and gauge structure of 
the low energy theory and another containing the tables relating the original 
lattice deformations with resulting gauge fields. 

2 Preliminaries 

We consider the tight-binding model in which the propagating electron dy- 
namics is described by the two-dimensional hexagonal lattice Hubbard model, 
given by the Hamiltonian, 

H=-tJ2 {aibn+a + Ca«n) , (2-1) 

where a}^ and the creation/annihilation operators for the electron on 

the A lattice site, and, respectively 6jj and 6„ are those for the B lattice sites. 
Each A site is connected to three nearest neighbor B sites by lattice vectors 
a = 1,2,3 (see the Figil]), 

i = a(l,0), 2 = a(-l/2,y3/2), 3 = a(-l/2, -73/2). (2.2) 

and vice versa each B site is connected to nearest A sites by lattice vectors 
— d. 

The lattice can be represented as a superposition of two simple (Bravais) 
lattices (A and B sublattices), each generated by lattice vectors {/i = 1 — 
2J2 = 1-3} (see the Fig. [I]), 

/i = a(|,-f), i, = a{lf). (2.3) 

The B-sublattice is shifted by 1 with the respect to the A-sublattice. The 
spacing for the Bravais lattice is |/i| = I/2I = y/Sa. 

The sublattice degeneracy has an important implication. Thus, a lattice 
function /(xn), will be described by a two component object fa,n, where 
a = A,B parameterizes the sublattice and n, Bravais lattice site within the 
sublattice. 

In the case of Hamiltonian fl2.ip . the electron's annihilation operators can 
also be combined into a two component field "^a,!! = (^n, &n)- In terms of the 
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Figure 1: The hexagonal lattice. The lattice vectors as seen from an A-site: 
At the ends of vectors a = 1, 2, 3 there is a B-site, and at the end of Bravais 
vectors k, i = 1,2 there is another A-site. 

new field the Hamiltonian ( 12. ip takes the following form, 

n 

^-tJ^K-D-'^n, (2.4) 



where 



D 



1 + r/ + r/ 





(2.5) 



where T^ , i = 1,2 represent elementary translation operators on the Bravais 
lattice. 

Taking into account the anticommutation relations, in terms of the new 
fields, 

{^a,n, ^Ji,^} = SapSn,m, (2.6) 
we can write down the discrete Lagrangian corresponding to the Hamiltonian 



eq. 



^ = iE<n■^.n-i^(^^^) 

an 

n 

+ t ■ ^B,n + ■ ^B,n-f, + " *B,„-[, + h-c)! (2.7) 
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The Low energy theory 

To find the low energy theory, let us go to the momentum space representa- 
tion of the action (\2.7\i . by using the Fourier transforms of the fields, 

$„(k) = ^vl>^,„e'''" (2.8) 

n 

where k = {k^,ky) is the lattice momentum. This definition implies, that 
the Fourier transform \l/a(k) is periodic with periods given by the reciprocal 
lattice vectors {ki,k2} satisfying, 

ki ■ ij = 27i5ij. (2.9) 

In our setup the reciprocal lattice vectors are given by 

h = ^{l,-V3), A;2 = 1^(1,73), (2.10) 

and = \k2\ = 4:71 /3a. 

The inverse Fourier transform is given by, 

'^a,n=^ [ d2fcv|/„(k)e-''^-", (2.11) 

^ JFD 

where the integration is done over the fundamental domain (Brillouin zone) 
-1/2 < ki < 1/2, i = 1,2 and Abz = Stt'^/^VSo^ is the area of the Brillouin 
zone. 

The action in terms of Fourier transforms takes the form, 

S=^[ dtd^k\i¥{k)^{k)+t¥{k)-D{k)-^{k)]. (2.12) 



where 



is the fermionic operator, and 

rf(A;) = 1 + e'*^-'"! + e*-'"^ (2.14) 

In the Brillouin zone the function d{k) vanishes for two distinguished 
values of k: = -K = (0,-^) and K2 = K = (0,^). These 
two points are called Dirac points. Since the Hamiltonian vanishes at these 
points, their neighborhoods are relevant for the low energy regime of the 
model. In the hexagonal Brillouin zone picture the Dirac points appear at 
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the corners of the zones. In the the rhombic Brillouin zone description, used 
by us, they appear in the core of the zone. 

As we restrict ourselves to the dynamics near the Dirac points, the leading 
contribution comes from the linear terms in the expansion of the fermionic 
operator, which is controlled by the expansion of the function (i(k), 

d{±K + k) = ± K-^ky)) + • • ■ = k ■ V(±K) + . . . , (2.15) 

where V(±k) is the momentum space gradient of d{k.). Using this expansion 
we can express the fermionic matrix D{k) in the following form, 

Di±K + A;) = f [i-it,)a2 T {-^ky)a,] + ... (2.16) 

The ± sign labeling the Dirac point. If we introduce an additional index 
labeling the Dirac point, then we can rewrite the fermionic operator in the 
leading order in the following form, 

D{k) = ^[{-ik,)f3 + {-iky)p] (2.17) 

where we introduced the matrices 

/3 = (T2®I, p = -(Ti®(T3. (2.18) 

If we chose 7" matrix as, 

7O = iag (g) as. (2.19) 
and introduce the other two Dirac gamma matrices as, 

= -7O/3 = -ai ® (T3, 7' = -7V = -^2 ® I, (2.20) 
as well as the Dirac conjugate of the wave function, 

^ = ^t.^o^ ^2.21) 

the action fl2.12p in low energy rewrites as, 

S = ^ [ dk^dky [i^(A;)7o^(A;) + ^i^{k){-ikj^^^{k)\ , (2.22) 

where the index m = x,y, and integration is performed in the vicinity of the 
Dirac points assuming some implicit energy cut-off. 

Representing now the low energy field as a Fourier transform of a contin- 
uous function \E'(a;), 



^(k) = I dh^{r)(^''-\ (2.23) 

2n 



where k is treated as an ordinary continuum space momentum, the action 
becomes, 



where i = 1,2 and vp = 3ta/2 is the Fermi velocity, which plays here the 
role of the speed of light. The Dirac Fermi field carries one index a = 
A, B labeling sublattice, one index ± denoting the Dirac point, in total four 
components. Here one can include also the electron spin index, which should 
rise the number of components up to eight, but we are not considering the 
spin component in this work. 

3 Defects and Hubbard model deformations 

So far we considered the two-dimensional lattice as a regular rigid object. In 
the real world, however, the lattice is a subject to smooth deformations as 
well as to discontinuous topological defects and impurities. 

The smooth deformations correspond to phonon fields. In the standard 
approach, the basic variable describing the phonon field is the position of the 
individual atoms in the lattice (see e.g. [IS!)- Then, in the linear response 
theory the phonon dynamics can be described by using the dynamical matrix 
encoding the elastic forces between the atoms. In the case of graphene is 
well established that the dynamical matrix should involve elastic forces for 
at least up to next-to-neighbor atom site, due to relative small difference in 
the distances to nearest and next-to-nearest neighbors. 

In the framework of the microscopic tight-binding model, the coupling to 
the fermionic field is realized through the point dependent modification of 
the hopping amplitudes in the Hubbard model. For small deformations the 
amplitude modifications are linear in the displacements of atoms from the 
equilibrium position. In this approximation, the description of the phonon 
field in terms of the displacements is equivalent to the description in terms of 
transition amplitude variations. For large deformations this may not remain 
true, however, the description in terms of amplitudes is advantaged because 
the amplitudes are always linearly coupled to the electronic field. In the low 
energy theory such deformations were shown to give rise to an interaction 
potential similar to gauge coupling [H]. 

On the other hand the study of various topological defects was shown to 
lead to the interaction potentials having the form of coupling to nonabelian 
gauge field backgrounds [T5HT9] . The distinguishing feature of a topological 
defect is that, in contrast to the phonon field, the 'defects field' may not 
necessarily have a limit as a continuous function in the low energy theory. 




(2.24) 
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This property of topological defects plays a crucial role in building up the 
nonahelian gauge field, since for coupling electronic modes at different Dirac 
points can be only achieved through discontinuous modes. 

All in all, both phonon fields and lattice defects can be taken into account 
in the tight-binding approach through the deformation of the Hubbard model 
Hamiltonian by adding to it certain types of operators. In what follows we 
will consider such deformations involving the local, nearest neighbor and 
next-to-nearest neighbor terms. Strictly speaking, certain types of topolog- 
ical defects should involve non-local terms of arbitrary long range but, as 
we can anticipate, the terms only up to the next-to leading are enough to 
generate a generic gauge field. 

3.1 Hubbard model deformations 

Let us turn to the tight-binding model described by eq. fl2.ip . and write down 
possible deformation terms. 

Consider first the local contribution. A generic lattice deformation may 
lead to local modification of the Fermi level. This modification can be taken 
into account by the following terms added to the Hamiltonian, 

Aifn = - ^(Qn^AnQn + h\^ZBj)n), (3.1) 
n 

with Zdni a = A,B are the components of a real lattice function. 

The nearest neighbor transition amplitude can also be modified. In gen- 
eral the amplitude can be complex, however, the corresponding term in the 
Hamiltonian should be Hermitian. Therefore, in the most general form it 
reads, 

Aifnn = - ^ (aiZn,abn+a + &|,+a^n,aan) , (3.2) 
n,a 

the bar stands for the complex conjugate. 

In addition to the nearest hopping the deformation, one can induce the 
next-to-nearest one. While for the undeformed model the amplitude of such 
hopping is expected to be much smaller than the nearest neighbor coun- 
terpart, its variation due to the lattice deformation could be rather large. 
This deformation is taken into account by adding the following term to the 
Hamiltonian, 

AH^nn = -Yl {(^i^An,ab(^n+a-b + ^La^ Bn,abK~b) ' (3-3) 
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where z^ah ~ ^nia^ a = A, B are the transition amphtudes for the hopping 
between nearest A-sites and nearest B-sites respectively. 

Thus, in total we have three types of parameters describing the lattice 
deformation: real valued site based (scalar) field Zn, link based (vector) field 
Zn,a and face based (2D-pseudoscalar) field ^n^g- 

Tout ensemble, the above modifications lead to the following deformed 
Lagrangian fl2.7p . 

L = Lq — Aifn — Aifnn " Ai/nnn 
n 

+ (^kn^i?n,i2^B,n-f, + ^kn^Bn,i3^B,n-f, + ^kn^Bn,23^B,n+;Vf, + h-c)] , 

(3.4) 

or in the compact form, 

n n,i=l,2 
n 

where we use the capitalized symbols to denote the following 2x2 matrices, 

Zrv = |(^A,n + ZB,n)^ + |(^A,n " ZB,n)(T3 + Z^^i^^CTi + 1(72) + Z^^i^{(Ti - 1(72) 

+ zl = zl zl = zl (3.6) 

Z^M = ^n.f.^+ + €j} + <,?.^3, = -z^ja- + -zlu^ + -zlt^s, (3.7) 

= + <[.^^3, = -zljjl + -zlu.a,, (3.8) 



with 



^n,h ~ ^n,2' ^n,i2 ~ ^n,3' ^n,h ~ ^",2' ^n,[2 ~ ^n,3' (3-9) 



2; 



2(^An,i2 + %n,i2)' /o ~ 2 (^An,i3 + ^Bn.is) ' (3.10) 



Z' 



^niii2 ~ i(^An,23 + %n23)' ^n,/i/2 ~ |(^An,23 ~ ^_Bn,23)- (3-12) 
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The Pauli matrices (Xi 2 and (J3, as well as unity matrix I are defined in 
the sublattice space. 

The 'lattice vector' Z^^, appearing in the Lagrangian fl3.5p . can be ex- 
pressed as a scalar product of a local vector field Zn and lattice vector /j: 

^„J. = Zn-l (3.13) 
The inverse transformation reads, 

Zn = 2^ E ^n,fi., (3.14) 
i 

where ki, i = 1,2 are the vectors of the dual basis given by eg. (12.91) . 
3.2 The low energy limit 

Let us now take the low energy limit of the theory described by the La- 
grangian (13. 5 p and find the coupling of the low energy electronic modes the 
modes of fields introduced in the previous section. In order to do this, let us 
again expand the action near the Dirac points for the fermionic field and keep 
only those modes of the deformation fields which couple to the low energy 
modes of the fermion. In the momentum space the expanded action reads, 

5 = 5o + ^ J d^kd''q[<i/l{k)Z{k - q)^+iq) + ^1(A;)Z(A; - q)^4q) 

+ ^l{k)Z4k - q)^-{q) + ^l{k)Z+{k - q)^+{q)] 
-i[^iik){Z{k - g) ■ V+ + Z*ik - q) ■ V_}vl/+(g) 
+ ^l(A;){Z(fc - g) ■ V_ + Z*{k - q) ■ V+}^_(g) 
+ ^\{k){Z^{k - g) • V- + Z*_{k - q) ■ V-]^-{q) 
+ ^>l{k){Z+{k - g) ■ V+ + ZX{k - g) • V+}M/+(g)] 
+ ^>l{k){e^-"-'Zi^^^{k - g) + e-^-"'Zl^,^{k - q)]^^{q) 
+ ¥_{k){e-'-'I^Zij^{k - g) + e'-"'Zl^^{k - q)]^ ^{q) 
+ ml{k){Z^,j,^_{k - g) + Zl^,^_{k - g)}vl/4g) 

+ ml{k){Zij^^{k - g) + Z^j^^ik - g)}M/+(g), (3.15) 

where the subscripts ± of Z±, Z± and Zi_^ indicate, that the field is eval- 
uated near the respective Dirac point, e.g., 

Z±(k) = Z(±K + k), (3.16) 
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absence of such a subscript means that the field is taken near the origin of 
momentum space. We also included the multiplicative coefficient into the 
momentum measure. 

In the matrix form the Lagrangian fl37[5|) takes the formQ 

L = Lo + Z^'^^^ -ai^cTj-^, (3.17) 

where /, J = 0, 1, 2, 3 and ctq = I, while the other sigma matrices are Pauli 
matrices. The field Z^'^ consists of three main contributions: 

Z^-^ = Z^-^ + Zll + ^nnn, (3.18) 

where Z^"^ comes from the local deformation of the Hamiltonian, Z^^ from 
the nearest neighbor and Z^^^ from the next to nearest neighbor deformations 
described respectively by terms AH^, Ai^nn and AH^ in the Hamiltonian. 
The tables of the values of ^nnnnnn terms of the original deformation 
parameters are given in the appendix [Bl 

Replacing the Hermitian conjugate wave function by the Dirac conjugate 
according to 'I' = \E'"'"7°, the Lagrangian fl3.18p becomes, 

L = Lo + + W^Ta^ + A^^-ff""^ + B^^^Y-Ta"^, (3-19) 

where the field $ and U = U"'Ta are Abelian scalar and su(2) algebra- 
valued scalar fields, while and = B'^Ta are, respectively, Abelian gauge 
(pseudo-electromegnetic) field and 5^ su(2) non-abelian gauge fields. The 
effective fields are expressed in terms of the original deformation modes as 
follows, 

* = -!(4-^1) + f(^"-^"')' (3-20a) 
= -\{zi_ + %_ - - %^) 

~ + ^x- ~ Zx+ — — i{Zy+ + Zy^ + + Zy_)], (3.20b) 

f^' = i(^i- + %- + ^i+ + %+) 

~ 8^[^a::- ~^ ~^ + Zx+ + i(^?/+ + Zy^ — Zy_ — Zy_)], (3.20c) 

= -\{ZA - zb) + ^iizl + zl) + {{z'^ + z''), (3.20d) 

for the scalar fields, 

A, = -\{ZA + zb) + l{zl + zl) - i(4° + 4°), (3.20e) 
Ai = ^(;2i-%)-f(z,.-zJ, (3.20f) 
A, = -^^{zy-Zy), (3.20g) 



^We are using the first factor for tlie sublattice space, while the second factor corre- 
sponds to the Dirac point space. 
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for the pseudo-electromagnetic field and, finally, 



+ 8^[^x'- ^ + Zx^ — Zx^ + 1(2;?/+ ~ Zyj^ — Zy^ + (3.20h) 
B\ = lizA+ + Za- + zb+ + Zb-) 

^jl^x- ~^ Zx- + Zxj^ + Zxj^ + "I" ~ ^?/- ~ ^J/-)] 

+ + + ' (3-20i) 

= \{za- — Zb- — za+ + zb+) 



2 — 4 

t[^X- + Zx- ~ ZxJ^ — Z^_^ — i{Zy_^ + Zy^ + Zy_ ^ ^J/ - ) ] 



+ l(/3+ ,3_^^3_^.3)^ (3.20j) 



Bq — 4(^1- ^i- ^i- 



^x- ~ Zx- — 2x+ + ~ KZy+ ~ Zyj^ + Zy- — (3.20k) 



^1 = ii.zA- - 2;a+ + Zb- - zb+) 

+ + (3.201) 

B\ = —\{za- — Zb- + za+ — zb+) 

"I" l\-^x- "I" ^x- "I" Zx+ + ^a;+ "I" "I" ^1/+ ~ ^y- ~ ^J/-)] 

+ i(z'_3 + 2:j^ + ^'3 + ^'3). (3.20m) 



i?o = !(S°-4)-¥(^'°-^'°)' (3-20n) 
Bl = l{zy + Zy), (3.2O0) 
^2 = + %) + K^x + Zx), (3.20p) 

for the nonabelian gauge field. The shorthand notations are explained in the 
Appendix [Bl 



4 Example: The St one- Wales defect 

As an example let us consider the Stone- Wales defect [20]. This defect con- 
sists of a couple of collateral heptagons paired by a couple of pentagons. As 
pentagons and heptagons are carrying opposite curvature 'charges', this de- 
fect appears like a two-dimensional gravitational vortex. As was found by 
microscopic simulations, this type of vortices are likely to be formed under 
strain as it relieves the stress in the direction of heptagons [21]. The defect 
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Figure 2: The Stone- Wales topological defect as a flip of two bonds adja- 
cent to our unit cell bond. The flipped bonds are of the same orientation. 
Choosing another pair of flipped bonds results in an equivalent configuration 
with A and B atoms in the turning bond exchanging the roles. There are 
three different possible orientations for Stone- Wales defect, corresponding to 
affected bond 1, 2 or 3. 



can be obtained from a regular hexagonal lattice by flipping two links adja- 
cent to a chosen (central) bond, as shown in the Figl2l As there are three 
inequivalent types of links connecting the neighbor A and B sites, there are 
three different orientations of the defect. Let us denote each orientation as 
type 1, 2 and 3 according to the orientation of the central bond of the defect. 

The Figl2] shows only the defect obtained from the alteration of a bond 
parallel to the unit cell, the other types are obtained through rotations by 
±120° of the shown one. Although, the differently oriented defects differ by 
just rotation of the hexagonal lattice, in the compact notations they have 
quite a different form. Apart from the small smooth part, the defect field 
can be represented in the Hamiltonian by nearest neighbor (nn) values +t 
corresponding to the two canceled dotted bonds and two next to nearest 
(nnn) values —t' ^ — t corresponding to the new emergent bonds. 

For definitiveness let us consider only the type 1 defect, drawn in the 
FigJ2J From the picture it is not difficult to deduce that the defect field is 
represented by the following (sub)lattice deformation field, 

where no is the location of the turned bond, which place the role of the center 
of the defect. Using the low energy expansion worked out in the previous 
section we find that, as a result of this defect the electronic wave function is 
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coupled to the following nonzero fields: The Yukawa scalar field, 
^ = ig[- cos(K ■ no + 7r/6)ri - sin(K ■ no + 7r/6)r2 + frg] ^(^^(r 

the Abelian gauge field given by, 

= (7(1/2, 0,0)5(^)(r- no), 
as well as nonabelian gauge field, 

= g{0, sin(K ■ no + n/6)n - cos(K ■ no + 7r/6)r2 + {V^/2)t^, 
— -\/3sin(K • no + 7r/6)ri + V3cos(K • no + 7r/6)r2 + {1/2)t^ 

x5(2)(r-no), (4.2c) 

where g = 3\/3ta^ /8n'^ = ^vpa. 

As can be deduced from the value of the parameter g, a single defect gen- 
erates vanishingly weak effective field. Therefore, a sizable effect is reached, 
when a finite density of Stone- Wales defects is considered. 

On the other hand, as the explicit form Eg. (14. 2 p shows up, the gauge 
and scalar depend discontinuously on the position of the defect, even for 
a fixed orientation. Indeed, since K ■ no = ^{tlqi — no2), where noi, for 
i = 1,2 are the integer lattice numbers of the defect position, there are three 
distinct values for the gauge and Yukawa field for this orientation of the 
defect, depending on (noi — no2) = 0, 1, 2 mod 3. A random distribution of 
Stone- Wales defects will result in the fields $ and 5^ with vanishing average 
Ti and T2 components, but non-trivial Abelian gauge field and components. 
Different results one obtains when the defects form a sublattice with fixed 
9 = K-n. 

5 Conclusion 

In this work we considered deformations of the hexagonal lattice Hubbard 
model describing the graphene's electronic wave function in the tight-binding 
approximation. These deformations involved arbitrary space and time depen- 
dent operators including local, nearest neighbor and next-to-nearest neigh- 
bor interactions. The low energy effect of the deformation is the Yukawa and 
gauge couplings to the electronic Dirac field. The low energy Lagrangian is 
given by the eq. (I3.19P and the fields are listed in eq. (I3.20p . There is no 
gravity-type coupling in this model since there is no spin connection in 2 -f- 1 
dimensional Clifford algebra. 



- no), 
(4.2a) 

(4.2b) 
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The considered deformation operators do not exhaust all possible defor- 
mations in general, since one may include arbitrarily long-range interactions, 
however, they are enough to saturate generic gauge coupling: It is easy to 
check, e.g. by an inspection of the Tables [3l H] and [5l that excluding one of 
these contributions would result in an incomplete gauge field. On the other 
hand including longer range interaction will not add any new terms to the 
low energy action, but merely redefine (renormalize) the existent ones. Al- 
though the way the gauge fields are redefined is an interesting topic to study, 
in particular because there are topological defects which involve long range 
deformations. 

The emergence of the whole gauge coupling is restoring the pseudo-Lorenz 
and internal symmetry. Moreover, it is upgrading the global SU(2) internal 
symmetry to a local gauge invariance, provided the defect generated fields 
are prescribed the proper transformation rules. Of course, the relevance of 
such symmetry depends on the dynamics of the defect fields not studied here. 
Therefore the inclusion of the defect dynamics into the consideration is an- 
other interesting direction of the development of this work. In particular, 
it would be interesting to investigate the possibility of having a controlled 
vacuum expectation value for the scalar sector of the model. This would be 
interesting from the theoretical point of view, as well as open new perspec- 
tives for practical applications. 
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A Dirac and Pauli matrices 

Reducible Dirac matrices 

We use the following definitions for the Dirac matrices, 

7° = ids (g) 0-3, 7^ = -o"! (g) 0-3, 7^ = -0-2® I. (A.l) 
We use the following definitions of Pauli matrices 

-Co)- 

The Dirac matrices satisfy the 2 + 1-dimensional Clifford algebra, 

{7m, lu} = '2Vtiu, V/iu = diag(-l, +1, +1). (A.3) 

The above Dirac matrices realize a reducible representation. By a unitary 
transformation we can go to a basis in which the Dirac matrices take the 
following form, 

7° = i(T3 ® I, Y = ai^ I, (A.4) 

where i = 1,2. In what follows, we show that this basis can be obtained by 
properly labeling the outer products of Pauli matrices multiplied by —7°. 



SU(2) algebra generators 

The above degeneracy of the Dirac matrices can be described in terms of 
additional SU(2) symmetry, generated by the following matrices, 

ri = (T2®(Xi, T2 = Cr2®(T2, T3 = I®(T3. (A.5) 

Direct inspection shows that matrices T^, a = 1,2,3 form an su(2) algebra. 
They also commute with the Dirac matrices 7^^, 

[Ta, 7l = 0, [Ta, n] = i2eabcrc. (A.6) 

In the diagonal basis of flA.4p the matrices reduce to outer products 
of identity and Puali matrices, = I TqI^ 



The diagonal basis 

The complete basis for the space of 4 x 4 matrices is formed by matrices 
Jfj.'Ta, Ta and 1, where, by slight abuse of notations, we mean irreducible 
7-matrices. We list the basis in the Table [H 

^To avoid the confusion, when speaking about the representation of (IA.4I) . we use the 
notations Ta for the Pauh matrices in order to underhne their su(2) nature. 
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n 


T2 




I 


1(8)1 




£72 <8 (72 


I<8)0-3 


70 


ic73 CT3 


iCTl (g) (T2 


— ifJl (Tl 


io-3®I 


71 


-CTi (g) (73 


0-3 8) (72 


-CT3 (g) fJl 


-0-1 (8 I 


72 


-a2 0l 


-I (8 0-1 


-I® 0-2 


-172 (8) £73 



Table 1: The basis for the space of 4 x 4 matrices 1 = I ® I, r^, 7^, and 7^rc 



In addition, we need to know the backward 'dictionary' to translate an 
arbitrary product of unity and Pauli matrices multiplied by —7°. These data 
we summarize in the Table [2l 






I 


0-1 


0-2 


0-3 


I 


-i(T3 <8) (73 = -7*^ 


<73 ® 0-2 = 7-*-^ 


-(73 (8 (71 = 7^r2 


-i(73 (81 = -7"t3 




0-2 8) 0-3 = -7^r3 


1(72 8) CT2 = ir2 


— i(72 (8 (7i = —in 


(72 (8 I = -7^ 


0-2 


-CTi (80-3 = 7-'- 


-ifJi (8 (72 = -7"ri 


i(7l (8 fTi = -7^x2 


— i7i (8 I = 7-^T3 


<73 


-il (8 (T3 = -ITS 


I (8 (72 = -7^r2 


— I (8 (7l = 7^Ti 


-il (8 I = -i 



Table 2: List of outer products of Pauli matrices multiplied by 70, —7° ■ {<7i 
aj), I, J = 0,1,2,3, in terms of gamma matrices and internal symmetry 
generators. The first factor changes by the row and the second factor by the 
column. 



B Field Tables 

The low energy coupling to the fermionic field obtained in the main part of 
the paper has the general form, 

~ ® aj)^, /, J = 0,1, 2, 3. (B.l) 

Here we are expressing the fields Z^'^ in terms of the original defect fields 
2^11; ^n,i and Zy^^. FoT the convenience, we are giving these expressions in the 
form of tables separately for each contribution: local (n), nearest-neighbor 
(nn) and next-nearest-neighbor (nnn), i.e. such that Z^-^ = Z^"^ + Zl^l + Z^i^^^, 
respectively in the tables [3], H] and [51 Subscripts "+" or "—"carried by the 
fields in the tables below denote that the momentum space modes near the 
Dirac point ±Kjl No subscript ± means that the modes near the origin 

■^In fact the argument should be ±(Ki — K2) = t2K, however the perioiiicity impHes 
that this point is equivalent to ±K. 
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are involved. In the real space this modes are given by the inverse Fourier 
transform, 

z,{r) = ^^j d2/cz,(k)e-'''" (B.2) 

The fields appearing in the Table |5] are related to the original defect fields 
(clarify (13. 6p ) by z = where ki are vectors of the dual basis. In 

coordinates, the vector z is represented by, 

= (^2 + 2:3)' Zy = -\^{zi-z^). (B.3) 

In addition, we use the following notations for the fields appearing in the 
tables, 

= i(^Ai2 + ^Bi2 + ^Ai3 + ^Bia)' (B-4) 
~ ^"^(^yli2 + ^312 ~ ^Ai3 ^ ^Bis)' (B-5) 
— i(^Ai2 ^ %i2 + ^Ai3 ^ ^is)' (B-6) 

+ ^M3), (B.7) 



2' 

^ _^ _^ ^ ^ y ^ ^ 

y ~ 2 V^A12 ^_B12 ^A13 



as well as 



l(^A23 + ^^23)' (B-8) 



^" = ite-%23)- (B.9) 





I 


0"! 


0"2 


0-3 


I 


i(z^ + Zb) 


j(2:a+ + Za- + 2:5+ + 


j(za- - + Zb- - zb+) 







{{zi + zi) 


K^i- + %- + + %+) 


+ ~ %+) 





0-2 


K^i -%) 




—j{zi_ — %_ — + %_|.) 





0"3 


2(^A — Zb) 


i(2;A_ - Zb- + - ^5+) 


^(zA- - ^B- - Za+ + 






Table 3: The local contribution Z^^ . The first index is counting rows while 
the second is counting columns. The fields are listed according to the defini- 
tion 03J]) 
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